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Abstract 

Mnltifractional Brownian motion is an extension of the well-known 
fractional Brownian motion where the Holder regularity is allowed to vary 
along the paths. In this paper, two kind of multi-parameter extensions 
of mBm are studied: one is isotropic while the other is not. For each of 
these processes, a moving average representation, a harmonizable repre¬ 
sentation, and the covariance structure are given. 

The Holder regularity is then studied. In particular, the case of an ir¬ 
regular exponent function H is investigated. In this situation, the almost 
sure pointwise and local Holder exponents of the multi-parameter mBm 
are proved to be equal to the correspondent exponents of H. Eventually, 
a local asymptotic self-similarity property is proved. The limit process 
can be another process than fBm. 

AMS classiEcation: 62 G 05, 60 G15, 60 G17, 60 G18. 

Keywords: fractional Brownian motion, Gaussian processes. Holder regular¬ 
ity, local asymptotic self-similarity, multi-parameter processes. 

1 Introduction 

In many applications, fractional Brownian motion (fBm) seems to fit very well 
to random phenomena. Recall that it can be defined by one of the four following 
properties. Let H G (0,1) [H is sometimes called the Hurst parameter). 

• is a centered Gaussian process such that 

Vs,t G R+; E [Hf i -\t- s\^^] 
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the process such that 

/■O 


Vt G R+; Rf = 


- {-u) 


H-i 


u)+ {t—u)^ ^.W{du) 

Jo 


is a fBm, 

• the process such that 

Vt G R+; Rf = [ 

Jr 1^1+2 

is a fBm, 

• B^ is the unique self-similar Gaussian process with stationary increments. 


Its efficiency has already been shown in simulation of traffic on Internet or 
in finance. This induced some recent progress such as stochastic integration 
against fBm. 

However, the main limitation of fBm is that the Holder regularity is constant 
along the paths. 

Multifractional Brownian motion (mBm) has been independently introduced in 
m and [ig. This process is a generalization of fractional Brownian motion 
where the Hurst parameter H is substituted by a function t ^ H{t). As a 
consequence the Holder exponent is allowed to vary along trajectories. 

The different definitions by the two groups of authors provided two different 
representations of mBm. 

Peltier and Levy-Vehel (ca) defined the mBm from the moving average 
definition of the fractional Brownian motion 

/•O _ _ /.i 








Jo 


where is a Holder function. 

Benassi, Jaffard and Roux (^) defined the mBm from the harmonizable 
representation of the fBm 




r - I 


.W(dO 


These two definitions were proved to be equivalent up to a multiplicative 
deterministic function ((HI). 

Moreover, in [g the covariance function of this Gaussian process has been 
proved to be 


E 


D{H{s),H{t)) 






where R is a known deterministic function. 

The goal of this paper is to study some multi-parameter extension of the mul¬ 
tifractional Brownian motion, ie a stochastic process indexed by R+, which is 
an mBm when = I. One extension has already been considered in 
2D extension of fractional Brownian motion has been already used in various 
applications such as underwater terrain modeling (El)- It may be more real¬ 
istic to allow local regularity to vary at each point : our extension of mBm in 
R^ may be used for this kind of application. 
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2 Multi-parameter extension of the fractional 
Brownian motion 

Since multifractional Brownian motion is an extension of fractional Brownian 
motion, we start with a review of the existing extensions of fBm. Most of the 
results in this section are well-known, but we give new proofs based only on the 
covariance functions. 

In the same way as Brownian motion has two main multi-parameter ex¬ 
tensions: Levy Brownian motion and Brownian sheet, two different multi¬ 
parameter extensions of fractional Brownian motion have been defined. 


2.1 Levy fractional Brownian motion 

This process can be seen as an isotropic extension of the fractional Brownian 
motion. Indeed, for the fBm, we have for all s,t G R+ 

E[Xt-X,f = \t-s\‘^^ 

A natural idea to extend this process for a set of index T C R+ is to 
substitute the absolute value by a norm. We get the Levy fractional Brownian 
motion, which is defined to be a centered Gaussian process of covariance function 

E[XsX,] = ^[\\sr^ +\\t\\^^ -\\t-sr^] ( 1 ) 

There are several definitions of this process by its trajectories. Among these, 
it can be defined as integral against white noise. Lindstrom stated the following 
(see 0). 

Proposition 1 The process defined by 


Xt = 



It - u\ 





W{du) 


is a Levy fractional Brownian motion up to a multiplicative constant. 


( 2 ) 


Proof This process is obviously Gaussian and centered. Thus we only have 
to show that the covariance function is of the form O- We have 


E[Xs-Xt]^ 



\t - u\ 





.du 


It — s — u\ 





.du 


We consider the change of variables from R^ into itself, v = where f is 

the linear application which maps the canonic basis of R^ to the orthonormal 
basis ^ei = pEfij: £ 2 , ■ ■ ■, ■ The differential of (f in any u € R^ is itself and 

the Jacobian 


•^<^(. 1 ......^) = ^)| =1 
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because the matrix of (j) is orthogonal. 
We have 


||t — s — u|p = (t — s — = ||t — s|p — 2 < t — S, M > +||u|| 

= ||t-sf-2||f-s||.ui + ||uf 


= (p-s|l-ei 


We obtain 


/R^ 


E[X,-X,Y = II 


.dv 


= 


\2H-N 




|ei 


t-s 


2 — 




P-5|| 


.dv 


and after the second change of variables, 

V = \\t— s||.w = p — s||/d.u> 


we get 


E [Xs - Xt]^ 







.dw 


therefore 


E = Km,h [Pf ^ + Pf" - P - sf"] 


□ 


The harmonizable representation of fractional Brownian motion can also be 
generalized. Before that, let’s recall briefly definitions of white noise and its 
Fourier transform. 

In the following, we will denote L^(R^) the set of functions / : > C such 

that \f{u)\^du < oo. 

Definition 1 The complex isonormal process is defined to be a centered Gaus¬ 
sian process W = |VF(/);/ G such that 


yf,g&Ll{R^)- E 


W{f)W{g) 


Then, white noise W can be defined by 



f{u)g{u).du 


W{E) = W{1e) 


Definition 2 A Gaussian process {wifyj & is said to be the Fourier 

transform of a complex isonormal process |vF(/);/ € L^{R^ ^)} if for all f € 

LUR^) 

Wif) = W{f) 


where f if the Fourier transform of the function f. 

The Fourier transform of white noise is defined in the same way. 
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This complex measure is usually used to define the harmonizable represen¬ 
tation of fractional Brownian motion 


B? = 


ojti _ I 


/R 


lei 


H+h 


M{dO 


that can be generalized in the following. 


Proposition 2 The process defined by 


Xt = 




e»(‘.S> - 1 


MidO 


where W is the Fourier transform of white noise in , 

is a Levy fractional Brownian motion up to a multiplicative constant. 


(3) 


Proof As will be done for multifractional Brownian held, the Fourier trans¬ 
form of the kernel of representation could be directly computed. But as this 
representation dehnes a real centered Gaussian process, it is enough to show 
that the covariance function has the form O- 

For all t € R'^, let’s denote by fi the function f i—> and consider the 

centered Gaussian process X = = W{ft)fi G R+|- 

First of all, let’s show that, almost surely, W{ft) G R. 

In fact, using W{ft) = W{ft), showing that /t G R is sufficient. 

Indeed, by E[Im{W{f))f = J^NiIm{f)f = 0, / G R imply W{f) G R 
almost surely, and 


Im{ft{x)) = 


/R^ 


/ 


Im (e 


— iKx^u^ 


ft{u)) .du 


sin < t — x,u > 


.du 


= 0 


by parity. 

The process X is therefore real and its covariance function is 


E 


E 


W{fs)W{ft) 

(e*<®’^> - 1) - 1) 


\2H-^N 


.d^ 


L 


IICII 

R^ 


ll^ll 


2H+N 


.df 


Then we have to consider 3 integrals of the form ^|g^ 2 ii+w -.df. 

As in proposition^] for t G R^ hxed, consider the change of variables from R^ 
into itself, u = where is the linear application which maps the canonic 
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basis of R'^ to the orthonormal basis ^ei = liln-, 62 ,..., eAr j, 
Then, we get 

r 1 _ g*<i,5> r X _ gj|hll-“i 

L- = in- 

After the second change of variables 


V = ||t|| .u = ||f||/d.u 
dv = 11^11“^.C?M 


we get 


f 1 — 


\\t\\^H+N r l_e»i 

Pir VlMP™ 


C'Ar,H>0 


Proceeding the same way for the 2 other integrals, we can conclude 

E[XA] = Cn,h [||sf ^ + ||<F - \\t - s|F] 


which shows that the process 
Brownian motion. □ 






is a Levy fractional 


2.2 Fractional Brownian sheet 

On the contrary to the Levy fractional Brownian motion, this process is not 
isotropic. In particular, we can have different Hurst parameters in each of the 
N directions. 

For the fBm, we have for all s, t G R+ 

E[X,Xt] = i + -\t-s\^^] 

As in the definition of Brownian sheet, another way to generalize fBm is to set 
the covariance equal to the tensor product of one dimensional covariances. Then, 
fractional Brownian sheet (fBs) is defined to be a centered Gaussian process of 
covariance function 

E =n ^ (4) 

i=l ^ 

As in the isotropic case, this process has two different representations by its 
trajectories. 


Proposition 3 The process defined by 


Xt = 


IR‘ 


N 

n 


\Hi-h 



W{du) 


is a fractional Brownian sheet, up to a multiplicative constant. 
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Remark 1 In Pontier/Leger introduced another moving average represen¬ 
tation of fractional Brownian sheet. 


Xt = 


N 

n - {-uf) 


sHi-k 


Proof This process is obviously Gaussian and centered. Thus, we only need 
to show that its covariance function has the expected form. We compute 


E [W,W] 



\U - Uif' 


2 



.du 





.dui 


We can see that the factor corresponding to each i, is the covariance of a fBm 
with Hurst parameter Hi (or a Levy fractional Brownian motion with = 1). 
Then we have 

N 

E[XsXt] = l[Ki,m - \u - 

i=l 

□ 


This process also has an harmonizable representation, using the Fourier 


transform of the white noise in as in the previous paragraph. 

Proposition 4 For all t = {tf), consider the function fit such that for all f = 

(CO, 

^ itmU _ 1 


Mu)= n 

m=l IsmI 


The process defined by 

Xt = W 


N 

{(ft) = n 


gjtmSm _ I ^ 

-W{df) 


is a fractional Brownian sheet, up to a multiplicative constant. 

Proof As in the previous proposition, let’s compute the covariance function 
of this process. 


E[XsXt] = 




N 

n 

m—l 

N 


n 


1 ) - 1 ) 




.df 


'R 


= n [N 


l^r, 


\2Hrr 


\2Hm+l 


■d^rt 


+ \tr 


\2Hrr 


\tm 


\2H„ 


using the same argument of the previous proposition. □ 
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Remark 2 The processes defined in propositions [3 and ^ are proved to have 
the same law. In fact, as a particular case of proposition^^ they are indistin¬ 
guishable. 


2.3 Stationarity of increments and self similarity 

Let us start by recalling the notion of increments in R;^. 

For a function / : [0, l]'^ —> R and h S R, one usually define the progressive 
difference in direction by 

A + hcz) - f{x) a x,xhci e [0,1]^ 

\ 0 either 

and for h G R^ and A = (ii,... ,ik), 

^h.Af — ^hi., .ii f ‘ ‘ ‘ ,ik f 

Despite the temptation to define the increments by Xt — Xg as in one dimension, 
it is better to set 

^Xs,t = ^t-s,{l....,N)^s 

re{0,l}^ 

If there exists i G {1,..., N} such that Si = ti, we have XXg^t = 0. Then, we 
consider 

7 = {j = l,...,Ar; Sij^ti} 

and 

Xt-s.lXs= Y 

^■£{0,1}#^ 


2.3.1 Isotropic case 

In the isotropic case, the following extension of fBm’s properties are well known 
(see ini). 

Proposition 5 Let X = G -R+| be a Levy fractional Brownian motion. 

We have the two following properties for all h G and a > 0 

Xt+h-Xh = Xt-Xo 

X n^X 

A-at — a A.t 

where means equality of finite dimensional distributions. 

Proof For all s and t in R(J^, we have 

E - Xh){Xt+h - Xh)] = I {e [Xg+h - Xhf + E [Xt+h - Xh]"^ - E [Xt+h - Xg+hf 

= E[XgXt] 


32 


For self-similarity, we compute 


E [XasXat] 


i (e + E [Xatf - E [Xat - Xas]") 

i (llosf ^ -h Ijotf ^ - ||at - asf 
a'^^E 

E [a^Xs a^Xt] 


□ 

Proposition [3 implies the stationarity of increments 

Proposition 6 The increments of Levy fractional Brownian are stationary, ie 
for all h€ R+ 

XXfi,t+h —^ XXo,t 


Proof We fix and write 

AXh,t+h = Y. i^lh^+ra.U - Xh) 

rG{0,l}^-{0} 

then in the development of E [AXh,s+hAXh,t+h\, we only have terms of the 
form 

E [{^[hi+riSi]^ - Xh) {X[hi+piti\, - Xh)] = E [X[r,s,].X[p.t.]J 

using the previous proposition. Therefore we have 

E [AXh,s+hAXh,t+h] = E [AXo,.AXo,t] 

□ 


2.3.2 Non-isotropic case 

In the non-isotropic case, the properties of self-similarity and stationarity of in¬ 
crements have been stated by Leger/Pontier (cf |H]). Here, we give another proof 
based on the covariance function rather than the moving average representation. 

Proposition 7 Let X = |At;t G R+'^ be a fractional Brownian sheet. We 
have the two following properties for all h G and a > 0 


AXh,t+h = AAo.t 


(^) 

Xat = A* 
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Proof We consider N independent fBm of Hurst parameter Hi, 

and the process Y = ^Yt't G R+| such that Yt = H^i ■ We can see easily 
that X and Y have the same covariance function. The same result follows for 
the increments ^AXh,t+h', t G R-+1 and ^AYh^t+h', t G R+ |- As a consequence, 
from 


N 


ay 


h,t-\-h 




re{0,l}^ 
N 


n[ 


yd) 


then we have 


E [AXh^s-\-h^^h,t-\-h] = E [AYh^s-\-hAYh^t-\-h] 


N 

n^[( 


y(*) _ _ Y^^^^ 

^Si-\-hi -^hi J ^hi J 


= i;[AXo..AXo,t] 

For self-similarity, we compute for all a > 0 
^ 1 


E [XasXat] = ^ - \ati - asi\^^^) 


2=1 


N ^ 


= a'^'^^^^E[X,Xt] 


= E 


X, Xf 


□ 


Therefore, we can conclude that both extensions of fBm satisfy the properties 
of self-similarity and stationarity of increments. 


3 The multifractional Brownian motion’s case 

Once again, we can consider two different kinds of multi-parameter exten¬ 
sion of mBm : isotropic and anisotropic extension. Note, first of all, that 
mBm already has a multi-parameter extension. Indeed, the formulation of Be- 
nassi/Jaffard/Roux in ^ was done for t G R^. We will see that it can be 
considered as an isotropic extension. 

3.1 Isotropic extension 

To define an isotropic extension of the mBm, the natural way is to substitute 
the constant H of the moving average representation of the Levy fractional 
Brownian motion, with a function. 


34 








Definition 3 Let H : > (0,1) be a measurable function. The process 

S -R+l such that 


X,= 




( 6 ) 


is called multifractional Brownian field. 


We will show that this process is the same as the process defined by Be- 
nassi/Jaffard/Roux. This result generalizes on the equivalence stated in the 
case iV = 1 in [^. 

Proposition 8 Let H : > (0,1) be a measurable function. The process 

defined by 

r gd*.?) _ 1 ^ 

Jr^ ||^||«(*)+- 

is indistinguishable, up to a multiplicative deterministic function, from the pro¬ 
cess defined by 0). This formulation is the harmonizable representation of the 
multifractional Brownian field. 

Proof First of all, let us compute the Fourier transform of the function ||.||“. 

(rii.ir,<^) = (ii.ir,^) 


/ pir ( 

/ e-^<'^’*>ip{w).dw] 

/r" ^ 

Jr" J 


we consider the change of variables 

X ^ R^ X R^ 

{w,t) ^ {w,X = cj}{t)) 

where (j) is the linear application which maps the canonic basis of R'^ to the 
orthonormal basis ^ei = n^, 62 ,..., We get 


inr.T) =11 

Jr" Jb 

^ / / 

Jr" Jb 


|agiAi||u)| 


ip{w).dw.d\ 
dw.du 




\N 


/R"jR" llw^lr llw^ll 

using the change of variables {w, A) 1 —> {w, u = ||ui|| A). Then we have 

1 


{nr,T) = 




/R 


N W\ 


I a-\-N 


ip{w).dw 


Thus, 




Ack 

||r(;||“+^ 
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We use this result to calculate the Fourier transform of \\t — .||“ — ||.||“. We 
will use the following property : if g{u) = f{u — a) then g = 

T[ii(-.r-ii.in(.) = [e-<«>-i]ii^ 

We deduce from this 

_ \\H{t)-§ _ II ||if(t)-f 


and Vt G R^, we have almost surely 


(u) — I 


<t,v> _ 2 




^ -2 






■«) = ^mt) / . ■w(dg) 


using the fact we saw previously that the second integral is almost surely real. 
Therefore, by an argument of continuity, the result follows. □ 


This process is obviously a centered Gaussian process. It is thus of interest 
to study its covariance function. The following proposition is an extension of 
the case fV = 1 stated in |H]. 

Proposition 9 Let | be a multifractional Brownian field. There 

exists a deterministic function : R ^ R such that the covariance function 
of X can be written 


E [X,W] 


Df^{H{s)^H{t)) ||s|| 


H{s)+H(t) 


( 8 ) 


Proof The easiest way to show this result is to use the harmonizable repre¬ 
sentation. By definition of W, we have 



- 1 ) 


1) 


.d^ 


This integral has already been calculated for a Levy fractional Brownian motion 
with a parameter H = Then we have 


E [W,W] 


f f 1 - e™' 

\\u\\H{>^)+H{t)+N 


H{s)+H{t) 


Di,iH{s)+Hit)) 




with D^^{x) = fj 






.du □ 


3.2 Non isotropic extension 

Another way to extend the multifractional Brownian motion for a set of index 
included in R;J^, is to copy the definition of the Brownian sheet. 
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Definition 4 Let H : —> (0,1)'^ he a measurable function. The process 

G -R+l such that 


Xt = 



\Hiit) I _ lHi{t) I 


where W is the white noise, is called multifractional Brownian sheet (mBs). 


As in the case of the isotropic extension, there also exists a harmonizable 
representation of the mBs. 

Proposition 10 Let H : (0,1)'^ be a measurable function. For all 

t = we consider the function (j)t such that for all ^ = (^i), 


N 


Mu)= n 

The process defined by 

Xt=W {4>t) = 


_ 1 

^ Ilf + ^ 

m=l 


N 

n 


_ I 


oN If |ffm(t) + 3 

U |<;m| 


--W{d0 


is indistinguishable, up to a multiplicative deterministic function, from the pro¬ 
cess defined previously. This formulation is the harmonizable representation of 
the multifractional Brownian sheet. 


Proof We have already seen that for each m G {1,..., N} 


T 

Moreover, we compute 

N 




r 


n 




\m—l 


Therefore 

/ N 


N 


n 


e*tm. _ 1 


(C 


r ^ 

/ n \\t^ 

JR'V L 


-X _ I 2 ; 

m •^m 1‘^m I 


/R^ 

N 


N 

n< 

m—l 






'^m \ 


.dx 


.dx 


1 ^™-- 




iU) 






T =W 


■ N 

If [ 1 ^" 


|Rm(t)-5 _ I 


A(t) 
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We use the same arguments as in proposition |H1 to conclude. □ 


The following proposition shows that the covariance structure of multifrac- 
tional Brownian sheet, is a generalization of the fBs’s one. 

Proposition 11 Let S -R+j" ^6 ® multifractional Brownian sheet. There 

exists a deterministic function D‘^ : —> R such that 

E {H{S) + H{t)) n + \t^\Hr.{s)+H^it) _ 1^^ _ + 

m—1 

(9) 


Proof As usually, we use the harmonizable representation of the process 

E[XsXt] = 


/r'^ 


N 

n 


n - 1) - 1) 

11 I/- 


m—1 


I Cm I 


ciSmim _ 


1) 




1) 


/R 




■d^n 


We remark that the factor corresponding to each m, is the covariance of a 
multifractional Brownian motion, with has already been calculated. Therefore 
we have 


N 


E [A,At] = n 


m—1 


|Rm(s) + Rmb) 


|Rm(s) + R™(t) 


|tm 




□ 


Remark 3 The form of the previous covariance function gives the idea to con¬ 
sider the process A = |yt; t € R+| defined from N independent multifractional 
Brownian motions A*^®) with parameter Hi by 


Yt = X 


( 1 ) 

t(i) 


..A 


(AT) 

t(") 


Although Y is not a Gaussian process, it is easily seen that it has the same 
covariance function as a multifractional Brownian sheet. This remark will he 
often used in the following. 


4 Regularity 

A lot of properties are known about the regularity of the trajectories of Brownian 
motion and fractional Brownian motion. As we will see, in the case of the multi¬ 
parameter extension of the mBm, we have to make some assumptions about the 
regularity of H before studying the continuity of trajectories. In the definitions 
of mBm (cf Q] and 0]), the function H is supposed to be Holder continuous. 
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4.1 Continuity of the two extensions 

We first recall the Kolmogorov’s criterion. 

Theorem 1 (Kolmogorov) Consider a process X = ^Xt;t € such that 

there exists C > 0, p > 0 and j > N 

Ws,tGR^, E[\Xt-X,\P]<C\\t-sr (10) 

Then, there exists a modification Y = € 12+| of X that is Holder con¬ 

tinuous of any order q € (0, 

As usually, the quantity K [[ATt — is studied for s, t S [a,b] where a ^ b 

and then, a patching argument is used to extend to s,t € R+. 


4.1.1 Isotropic case 

Lemma 1 For all rj and p, such that 0 < p < pL < 1, the multiplicative fac¬ 
tor of covariance function in 01, is positive and belongs to C°° {[r], p]). 
Moreover, its order n derivative is given by 


D 




1 - e* 


I x-\-N 


In" 


T.du 


( 11 ) 


Proof As the integral of a positive function, is positive. By an argument 
of uniform convergence of integrals m on [? 7 ,/i], is C°° {[r],pi]) and the 
derivatives are obtained by derivations of the integrand. □ 


Proposition 12 For all s,t G [a, 6], we have 

]-E[Xt-Xsf = ll[il(s)+il(t)] X 


1 


+ - 


d'^ip 

dx"^ 


{His)+Hity,\\s\\ 


d'^ip 

dx'^ 


+Oa,b [{H{t) - H{s)) 
where p{x,y) = D{x)y^. 


- s 


(iJ(s)+il(t);||t||) 

Oa,b (Hit) - H{s)f 


X {Hit) - His))^ 

( 12 ) 


Proof Using the covariance function of the multifractional Brownian field, we 
have 

^E [\X, - Xt\‘^] = D[2il(s)] ||sf^(«)-i?[il(s) + il(t)] 

+D[2H{t)] -D[H{s)+H{t)] p||^W+^(‘) 

+D[H{s)+H{t)] (13) 

We have to get a second order expansion of this expression. 

We introduce the function p defined by 

p{x,y) = D{x)y'^ 
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We can write 


^E[\X,-Xt\^] = ^{2H{s),\\s\\)-^{H{s)+H{t),\\s\\) 

+>f{2H{t), ||t||) - >p{H{s) +H{t), ||t||) 
+D[H{s) + H{t)] (14) 


We use the second order expansion 
^{2H{s),\\s\\)-^{H{s)+H{t),\\s\\) 


{His)-Hit))x^iH{s)+H{t),\\s\\ 


+Oa,biHis) - H{t)f 




An inversion of roles between s and t provides the expansion of 


^{2H{t), ||t||) - tf{H{s) +H{t), 


Then 1141) becomes 
1 


E[\X,-Xt\^] = (Hit)-His)) X 


dip 


^ (Hit)-His)) 

' 7^ - ^ 


__dx 


iHis)+Hit), 


dp 

dx 


iHis)+Hit),\\s\\ 

av 


iHis) + Hit), ||a||) + ^ iHis) + Hit), ||t||) 


ED [His) + Hit)] lit - + Oa,b iHit) - His))'' 


Since 


iHit)-His))x 


dp 


iHis)+Hit), 


dp 


_dx ’ dx 

is Oa,b [{Hit) — His)) (||t|| — ||s||)], the result follows. □ 


iHis) + Hit),\\s\\ 


Corollary 1 For all s,t G [a,b], we have 

^E[Xt-X,]'' = D[2Hit)]x\\t-sf^^*^ 

+ ^i2Hity,\\t\\)x(.Hit)-His))' 

+Oa,b iHit) - His))'' + Oa,b (p - (15) 

where pix,y) = Dix)y^. 

Proof Using the expansion of D [His) + Hit)] and 

\\t - s\\^is)+Hit) ^ p _ ^||2ff(t) _ p _ ^11 ^ 

we get 

i^[if(s)+i7(t)] X = i^[2i7(t)] X (16) 

+Oa,b (||t - + Oa,b iHit) - His))'' 
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Moreover as Hit) <1 for all t G [a, b ], we have e = I — H{t) > 0 and 


2 {H{t) - H{s)) (pll - ||s||) = 2 {H{t) - H{s)) (p|| - ||s||)^ x (||t|| - ||s||)^-^ 

< (pil-piir + dltll-piD^- 

that implies 

iH{t) - His)) (lldl - ||s||) = (Hit) - His))^ + {\\t - sf «(*)) (17) 

We conclude by (tT^ . (tT?)ll and itTTIl using first order expansion of in x and 

y. □ 

Using the continuity of D, D' and D”, we can state from the previous propo¬ 
sition 

Corollary 2 There exist positive constants K and L such that 

Vs,te [a, 6 ]; E[Xt-Xsf <K\\t-s\\'^^^^'^+L \H{t) - His)\^ (18) 

Corollary 3 Suppose H is (3-Holder continuous. There exists a constant M 
such that 

Vs,te[a,5]; U [X* - < M ||t - sf ( 19 ) 


4.1.2 Non-isotropic case 

Lemma 2 There exists positive constants K and L such that 
ys,tG[a,b]; E[\Xt-X,\^] < if ||t - sf ^ + L ||i7(t) - iJ(s)f (20) 

Proof By remark|21 we have 


E[Xs-Xt]^ = E 


= E 


' N N ' 

n^2>-n^.'s 

fnM.’.-vsinv 

. \2=1 i>l 


.{i) 

s(0 


/N-l \ N 

nA,<sui~>-nvg, 

V 2=1 / 2=1 


TT TT ^( 2 ) 

LL^s(^) ^i(l)^£C2) 

2>1 2>2 / 

2 


By the inequality of convexity (X) Oi)^ < 


E [Xs -Xt]^ <N{e 


r 


2 




r 

nv:i 

E 

■y(l) y(l)' 

\--- + E 

■ (AT) y(Ar)' 

1 

.i>l 







'N-1 


1 2 

E 

n V:i 

I 


. i=l 

i 


Since there exists a constant M = Xla^b such that 

2 


Vt € [a,6], Vi; E 


X 


(i) 

t(0 


< M 
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we get 


N 


E [X, - Xtf < ^ E X% - X% 


( 21 ) 


Using 


^t(>) 


E 

eu implies 


< + L, {H,{s) - H,{t)f ; Vz = 1,..., TV 


E[Xt-Xs\ < NM'^' 


N 


N 




\i=l 


\i=l 


□ 


Corollary 4 Suppose E[ is P-Holder continuous. There exists a positive con¬ 
stant M such that 

Vs,te[a,6 ]; U [X* - < M||l - sf (22) 

4.1.3 Existence of a continuous modification 

In both isotropic and anisotropic cases, under Holder regularity assumptions for 
iJ, we have an inequality 

E[Xt-Xs]^ <K\\t-sr 

But to use the Kolmogorov criterion, we need to have a > N. 

As the random variable Xt — Xg is Gaussian, we can write, for each integer n 

and choose n such that n.a > N. 

We conclude by a classical patching argument. For a and b, Kolmogorov’s 
theorem gives a continuous process 1"“’^ = t G [a, 6]|. Consider a' and b' 

such that [a, 6] C [a',b']. The processes and 1"“ coincide on [a, 6]. Thus 

ytG[a,b]; = = i 

and, by continuity 

= [a, 6 ]} = 1 

Then we can define a process Y on who coincides with on [a, b] and 
we can seen easily that this process is continuous. 
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4.2 Holder exponents 

The notion of Holder function is well known. It is interesting to consider a 
localized version of this notion. 

For the paths of a process X, one usually define two kinds of exponent (see 

n, El): 

• the pointwise Holder exponent 




= sup < a; limsup 
I p^o 


P° 


Xt - X. 


< oo 


the local Holder exponent 


d(to) = sup ■ 


a; limsup sup 

p—*-0 s,t£B{to,p) 


\Xt - Xs 


< oo 


We can see easily that for all to, we have 

a{to) < a{to) (23) 

A study of these exponents, in the case of ID mBm, is made in [2]. 

Remark 4 If H is fS-Holder continuous, then the local Holder exponent (3{t) of 
H at every point is not smaller than fS. 

Conversely, suppose that the local Holder exponent of H at every point of 
a compact [a, b] is positive. Then H is (3-Hdlder continuous on [a, b] with (3 = 
inft6[a,b]/3(i)- 

In the following, we suppose that H admits positive local Holder exponent 
ft (to) at every point t. 

Proposition 13 Let X = € -R+| be a multifractional Brownian field. 

For all to € R+, the local Holder exponent of X at to is almost surely given by 

a{to) = PHo) A H{to) (24) 

and the pointwise Holder exponent of X at to satisfies almost surely 

aito) = Pito) A H{to) (25) 

where P{to) and (3{to) denote the pointwise and local Holder exponents of H at 
to- 


As a consequence of this result, if H satisfies 

VteR)^; (3{t) <H{t) 

the Holder regularity of multifractional Brownian field of parameter function 
H is given by the regularity of H (and not by the value of H). This point is 
developed in [7]. 
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Proposition 14 Let X = | be. a multifractional Brownian sheet. 

For all to G -R+, the local Holder exponent of X at to is almost surely given by 

a(to) =/3(to) A miniJi(to) (26) 

i 

and the pointwise Holder exponent of X at to satisfies almost surely 

a{to) = P{to) AmmH^{to) (27) 

i 

where /3(to) eind (3{to) denote the pointwise and local Holder exponents of H at 
to- 

The proofs of propositions na and m are detailed in the three following para¬ 
graphs. 

4.2.1 Lower bound for the local Holder exponent 

A lower bound for the local Holder exponent is directly given by Kolmogorov’s 
theorem. Indeed, for X a multifractional Brownian held or a multifractional 
Brownian sheet indexed by [a, 6 ], for all n € N, there exists A„ > 0 such that 

with a = 2 inf[Q ;i](/3 A H) or a = 2 inff^ {,] (/3 A min^ Hi). 

Kolmogorov’s theorem states that there exists a modihcation of X, which is 
g-H61der continuous for all g G (Ojf ~ Then, for all to G R+ and all 
a, 6 G R+ such that a <b and to G (a, b), we have 

V / TVT ~ \ Ck N 

Vn G N; a{to) > - - — 

2 Zn 

and therefore, taking the limit n ^ oo 

a{to) > 

As H is continuous, we can take the limit (a, 6 ) —> {to, to) and we get 


• in the isotropic case. 


a{to) > P{to) A H{to) 

(28) 

• in the non-isotropic case, 


a(to) > /3(io) A minRj(to) 

(29) 


I 


4.2.2 Lower bound for the pointwise Holder exponent 

By Il2.‘lll . nara.granh 14.2.11 nrovid es a lower bound for the pointwise Holder expo¬ 
nent. However, it can be improved in the case /3(to) < Pito)- 
Let X = G R+| be a multifractional Brownian held. By corollary |21 

there exist positive constants K and L such that for all s,t G R(^, 

E [Xt - Xsf <K\\t- sf + L \H{t) - R(s)|2 
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and by corollary O there exists positive constants a and M such that 


[a, 6]; E[Xt - < M\\t - sr 

Therefore, using Kolmogorov’s criterion, there exists a modification of X, which 
is j^-H61der continuous for all i/ S ]0, y[- In the following, we consider such a v 
with i € N 

For all e > 0, there exist po > 0 and M > 0 such that for all p < po and all 
t € B{to,p) 


E 


Xf-X, 


1 2 


to 


<M p^ 


Then, setting 7 = /3(to) A (to) — C; for all p S N 


P{\Xt-Xt,\> p^}<E 
Let p = and for all to S N, 


Xt - X,, 


-|2p 


< Mp pP^ 


Dm = {to + fc.2-fo+’"); fc e {0, ±1,..., ±2™}"^} 


let us compute 


P 

< 


|-^to+fe.2-(™+") ^to\ 

max -—_— - > 1 

fcG{±l,...,± 2 ™}"‘ 


Z P{IX, 

fee{±i,...,± 2 ™}'^ 


2 ”7?T 


-x,j > 2 -^-} 


< Mp 2 "^'’” 


Let us take to = = nn and p € N such that — pe < 0. By the 

Borel-Cantelli lemma, there exists a finite random variable n* such that almost 
surely, 

Vn > n*; max iXj +fe. 2 -(i+'")" “ ^tol < (30) 

fce{o,....±2''"}'^ 

From (EOll, we show that, almost surely, for all m G N, we have 

ytGDm; \Xt-XtJ<C2-'^’- (31) 

• if 0 < TO < nn, ED follows directly from 

• if TO > KTi, for t G Dm, let 


Cp," = {a; G Di^ri',yi, (to)* < Xi < U} 

Then consider t € B{t, n Cj™*. 

As the paths of X are ^-Holder continuous, we have 

|X( -Xt\< C2-^(i+fo" < C'2-T'" 


and by (Tinil . 


|X,-- XtJ < 2 -^" 


Using the triangular inequality, the result follows. 
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Therefore, ED leads to 


yme'N-,ys,te D^; jX* - < 2 ^ 2 "^" 

Using the continuity of X and m +oo, we get 

sup \Xt - X,\ < 2C 

s.ieB(to, 2 -") 

and therefore, almost surely, 

r ^ , 

iimsup sup - < +00 (32) 

p—>0 s,teB{to,p) P~^ 

By (TT^ . for all e > 0, almost surely 

Q;(to) > l3{to) A i?(to) - e 

Taking e G Q+, we have almost surely 

a{to) > Pito)H{to) (33) 

For a multifractional Brownian sheet X, by lemma |21 we get in the same way 
that, almost surely 

a(to) > Pita) A i^i(^o) (34) 

for alH = 1 ,..., X. 

4.2.3 Upper bound for the pointwise Holder exponent 

The main result getting the upper bound for the Holder exponents, is the fol¬ 
lowing lemma, a direct consequence of proposition^] using continuity of U, D' 
and D". 

Lemma 3 Let X = t G -R+| be a multifractional Brownian field. For all 
[a, 6 ] C , there exist positive constants ki, k 2 , h, h such that 

Vs,tG[a,6]; E [Xt - X,f > k^ \\t - - h {H{t) - H{s)f (35) 

E [Xt - Xsf > k 2 {H{t) - H{s)f - h \\t - sf (36) 

Proof We only have to study the multiplicative factors of \\t — and 

{H{t) — -ff(s))^ in itT^ 

• Let fci = inftg[a_b] D [2H{t)] and h = suptg[„ ^ D [2H{t)]. 

By continuity oit^D [2iJ(t)] on the compact [a, 6 ] and as the function 
D is positive (lemma P), for all t G [a, 5] 

0 < ki < D [2H(t)] < h < +00 


• and let 

$(t) = ||tf X {D [2H{t)] In^ ||t|| - 2D' [2H{t)] In ||t|| + D" [2H{t)]} 
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By lemma n 


r 1 _ 

' -V-' 

>0 

Let k 2 = inftg[o,b] ^{t) and h = sup^gj^ ^ $(<). 

As previously, for all t S [a, 6], we have 

0 < k2 < ^{t) < h < +00 


□ 


Lemma 4 Let X = € -R+j ® multifractional Brownian sheet. For all 

[a, 6] C -R+, there exist positive constants ki, k 2 , li, h such that 

Vs,ts[a,6]; t — sGR+.Ci 

E [Xt - Xsf > fci lit - sf - h {m{t) - H,{s)f (37) 
E [Xt - X,f > k2 {H,{t) - H,is)f -hWt- sf (38) 


Proof For all s, t such that t — s G R+.e^, using lemma |3 we have 


E[Xt-Xs]^ = E 


yi'^) ylL 

^t (0 


U) 

tO) 


> h -hiH,{t)-H,is)f 


and 


E [Xt - Xsf > k2 {H,{t) - H,is)f - h [U - 


□ 


From this result, the upper bound for the pointwise exponent is a conse¬ 
quence of the following lemma whose proof is the same as the case A^ = 1 (see 

m) 

Lemma 5 Let X = G -R+^ be a Gaussian process. Assume there exists 

pi G (0,1) such that for all e > 0, there exist a sequence of ) 

converging to 0, and a constant c > 0 such that 

\JnGN-E [Xt+h^ - Xtf > c||h„f 

Then we have almost surely 

a{t) < p. 

Let X = ^Xpt G R+ I be a multifractional Brownian field (resp. multifrac¬ 
tional Brownian sheet). Let /3(to) be the pointwise Holder exponent of H at tg- 
We consider the two cases : 
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• if H{to) < /9(io) (resp. Hi{to) < l3{to)), by definition of /3(to), we have 
,■ \\H{to + h)- H{to)\\ _ 

Hence, by (resp. lITTll i. there exists a positive constant C such that 


Then, by lemma 0 


a(to) < H{to) (resp. Hi{to) ) (39) 

• if iJ(to) > /3(io) (resp.iJi(to) > /3(<o)), we consider a e {/3(to)] H{to)) 
(resp. a G Hi{to))). There exists a positive constant C and a 

sequence (/in)„g]\j converging to 0 such that 

Vn e N; \\H{to + K) - i?(to)|| > C\\Kr 

Then, by IpM (resp. l|2Hl)) 

E[Xt,+^-Xt,f > 


hence, by lemma 0 


a > a(to) 


and therefore 

a{to) < P{to) 


(40) 


We can restate the upper bounds (P|l and iBUI) of the pointwise Holder exponent 
of X at to 

a{to) < P(to) A H{to) (resp. /3(to) A H,{to) ) (41) 


4.2.4 Upper bound for the local Holder exponent 

By (I23II . any upper bound for the pointwise Holder exponent is an upper bound 
for the local Holder exponent. But we can improve on this result in the case 
P{to) < H{to). We first give an analogous of lemmaj^for the local exponent 

Lemma 6 Let X = ® Gaussian process. Assume there exists 

pL € (0,1) such that for all e > 0, there exist two sequences {hn).^^pf and iln)n£N 
of converging to 0, and a constant c > 0 such that 

WnGN; E [Xt„+h^ - Xt„+iJ^ > c\\h^ - 
Then we have almost surely 

a(to) < pL 
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Proof Let 6 0 and considGr two scqucncGS (^n)nGN and (^n)nGN Q'S in tliG 

statement. 

For all n G N, the law of the random variable is A/’(0, cr^). 

From the assumption, we have an —> +oo as n —>■ +oo. 

Then, for all A > 0, 


P 


\\hn-lnr+^ 



P 


l^tp+hr, — ^tp+Zn I i_ 

\\hn-lnV+'' A 



1 

^/2^1Tan 


exp 




.dx 


+ + 00 


1 


Therefore the sequence -Xt +; | ) converges to 0 in probability. 

then there exists a subsequence which converges to 0 almost surely. Then we 
have almost surely a(to) < A* + £• Taking e € Q^, the result follows. 

□ 


Let a G {Pita); H(to)) (resp. a G (/3(to); iLi(to)))- As 

\H{t)-H{s)\ 

hmsup sup —n-n-= +oo 

for all M > 0, there exists po > 0 such that 

Vp < po;3s,t G B{to,p); \H{t) - H{s)\ > M||t- s||“ 

Therefore we can construct two sequences (/i„) and (A„) converging to 0 such 
that 


Vn G N; |iL(to + ^n) — H{to + ln)\ > M\\hn — InW^ 

By lemma El we can deduce 

a{to) < /3(to) (42) 

4.3 Directional Holder exponents 

One may also define directional pointwise and local Holder exponents in the 
direction u G = |u G R^; ||m|| = l| by 

a„(to) = sup I a; lim = qI 

P^O p°‘ j 

and 

~ ^ f r \Xt-Xs\ ^ 

<au{toj = sup < a; hmsup sup -r-^— < oo > 

I p^o s,tGB{to,p) Ir ~ 'SIp 
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As previously, for all u G iY, we have 


au{to) < au{to) (43) 

Moreover, we can see easily that for all u GU, we have 

a{to) < auito) and a(to) < a„(to) (44) 

Proposition 15 Let X = ® multifractional Brownian field. 

For all to G and all u G U, the local Holder exponent of X at to in the 
direction u is almost surely given by 

auito) = Pu{to)/\ H(to) (45) 

and the pointwise Holder exponent of X at to in the direction u satisfies almost 
surely 

aufto) = Puito) ^ H(to) (46) 

where Puito) o-'^d Puito) denote the pointwise and local Holder exponents of H 
at to in the direction u. 


Proof Let to G R(^, u GU and consider the stochastic process 

^ = Xtg+p,u\P > 0 | 

By definition, auito) and ctuito) are respectively the pointwise and local Holder 
exponents of X at 0. 

Let Hip) = Hito + p.u). We have 


-E 

2 


Xp-Xr^ 




= D 


2iL(p)] \p - (2iL(p); ||to + P-^ll) x (iJ(p) 

+o(ff(p)-R(??)) + o {\p - 


by corollary n 

Then, using the same method as in proposition II61 the result follows. □ 


Proposition 16 Let X = G be a multifractional Brownian sheet. 

For all to G R^, the local Holder exponent of X at to in the direction Ci is 
almost surely given by 

a-ePto) = PePto) L. Hiito) (47) 

and the pointwise Holder exponent of X at to in the direction u satisfies almost 
surely 

UePto) = PePto) L. Hpto) (48) 

where P^pto) and P^pto) denote the pointwise and local Holder exponents of H 
at to in the direction e^. 
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Proof As in the proof of lemma |21 there exists a constant M > 0 and a 
one-parameter mBm such that 


= M.E 




*0 +P-<^i 


-X 


*0 +'?■<=» 


Then, using the same method as in oronosition II41 the result follows. □ 


4.4 Application of Dudley’s theory 

Another way to study the regularity of our processes is to examine the behavior 
around zero of the modulus of continuity 

(Ax,TiS) = sup \Xs-Xt\ 

s ,t^T ‘,d{s ,t) <S 

When the process studied is Gaussian, it is convenient to consider the 
pseudo-metric 

dis,t)=E[X,-Xt]^ 

As usually, we define the ball of radius r > 0 about t G T by 
Bd{t,r) = {s € T;d{s,t) < r} 

and we say that (T, d) is totally bounded if for all e > 0, there exists E,... ,tjn G 
T such that 

m 

TC \jBd{k,e) 

i=l 

When (T, d) is totally bounded, we can define the metric entropy e i—*■ D{e, T, d) 
where D(e, T, d) is minimum number of balls of radius e required to cover T. 

The following theorem allows to improve on the results of the previous para¬ 
graph. 


Theorem 2 (Dudley’s Theorem) Consider a centered Gaussian process X = 
G indexed by the pseudo-metric space {T^d). If (T,d) is totally 

bounded and if a/I n D{r, T, d)dr < oo, then X has a continuous modification 

Y = ■! Yt;t G 1 . Moreover, there exists a universal constant C > 0 such that 


lim sup 
( 5^0 + 


_ tAy.TjS) _ 

/o A/lnD(§,T, d)dr + C(50nlni 


< 24 


To apply this result, we first need to verify the assumptions about the metric 
entropy. 

Lemma 7 Let T C R^ measurable and d a pseudo-metric on T. If there exists 
C > 0 and a > 0 such that 


ys,teT-, d{s,t) < Clls-fir 
then there exists rp > 0 such that Vr G [0, rp], 

D{r,T,d) < C^.Leb{T).r-^ 
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We saw previously that the 2 multi-parameter extensions of the mBm, satisfy 
the assumption of this lemma with T = [a, 6]. Then there exists Ca^b > 0, 
a = aa,b > 0 and rg > 0 such that for all r G [0, tg], 

D{r, [a,b],d) < .Leb{[a,b]).r~'^ 


As a consequence, ([a, b],d) is totally bounded and in the neighborhood of 0, we 
have _ 


Win £)(r, T, d) < \ K - In? 

V a 


therefore the integral ^yhiU(r^T\^dr is finite and we can apply Dudley’s 
theorem. We get 


Proposition 17 Let X = G -R+| be one of the multi-parameter exten¬ 

sion of the mBm. For all a ^ b, there exists Ca.b > 0 and a = aa.b > 0 such 


that 


lim sup 
e^O 


sup||s-t||<, 

hW 




< Ca,bVN 


Proof First of all, we study the quotient 

/o v'lnD(§,r, d)dr + CJ^lnlni 


N{S) 

m 


The derivative of the numerator is 


NfS) = ^\nDi^,T,d) + C y In Ini 


21ni^lnlni_ 


< \l K -In (5 -|- C\/ In In 


1 


1 - 


2 In I In In j- 


and the derivative of the denominator is 


D'{S) = Win 


1 - 


21ni 


Then we have 


V N'{S) 

- 

5^0+ D'{d) 


and by a L’Hopital’s rule type argument. 


lim sup 
( 5 ^ 0 + 


N{S) 

m 


< 



52 



























Then we have 


sup^(s |Xs - X(| 
lim sup- -== - 

The problem is now to transform sup^^^ into sup||s_(||<^. 
To do this, we write 



Vs,t € [a,b]-, d{s,t) < Ca,b\\s - t\\ 


then ||s — t|| < e = implies d{s, t) < 6 and we get 


SUP||^_t||<JXs - Xjl 

hm sup-- ' -< 24 

Ca,be^^a\n\ 

which gives the expected result. □ 


a 



This result is more powerful than knowledge of Holder exponents. It gives 
the behavior oi \Xs — Xt \ in a ball around to- 


5 Locally asymptotic self-similarity 


Extending fBm into multifractional Brownian motion implies the loss of the two 
properties of self-similarity and stationarity of increments. However, a weak 
form of self-similarity remains, called locally asymptotic self-similarity (see pQ, 
^). As we will see, this property still holds for the two kinds of extension of 
mBm in R'^. 

Theorem 3 Let X = G -K+^ ^6 o multifractional Brownian field. 

For allto € the law of the process (p) = |yf(p) = ^*o+p" 
converges weakly if one of the following two conditions holds 

1. a = H{to) and H{to) < infu_.„ fiuvitfij 

where fiuvito) = sup |a; limp^o = o|. 

Then, the limit measure is the law of a fractional Brownian field with 
parameter Hitfi). 

2. a = vaiu,v (iuv{to), H{to) > iniu,v Puv{to) and for all u,v G R+, the 
following limit exists 


\H{to + pu) - g(to -k pv)\ 
p™ pinf„,„/3„„(to) 


T{u,v) 


with (u, v) 1 -^ bounded on [a, 6]^ for some (3 > 0. 

The limit measure is the law of a Gaussian process guch ikat 


E 


Y™G,v duv{to) _ Y™G,v huv(to) 

^ U ^ V 


2 


Kto [r(M,n)]^ 
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Remark 5 As in the Levy fBm’s case in proposition\^ the same result as 
theorem\^can be stated for the increments AX defined in section m The law 
of the process Y°‘{p) = ; u G -R+| converges weakly under 

the same assumptions. 


In the case = 1, for all u,v G R+, we have f3uv{to) = P{to)- Therefore, 
theorem 0 has a simpler statement. The two cases to be considered, depend of 
the comparison between R(to) £md the pointwise exponent /3(to) of H. 

The following example shows that the limit considered in the second case, 
can be non trivial. 

Example 1 In the case = 1, let H{t) = j + t^ for t G [0, j]. 

For to = 0, we compute, for all u, v and p > 0 


\H{p.u) - H{p.v)\ 
1 

P2 


1 1 




< 


\u — n| 


1. 

2 


The limit measure is the law of a centered Gaussian process Y such that 
E[Yu-Y„f = Ko 


ie 

E[YuY.,]=Kou^v"^ 


Theorem 4 Let X = G R+| be a multifractional Brownian sheet. 

The law of the process Y°‘(p) = i u G R+| converges weakly 

if for all i G {!,..., A^}, one of the following two conditions holds 


1. ai = Hi{to) and Hifto) < in.iu,v Plyito) 

where Puy{to) = sup |a; limp^o \^dto+pu)-jii{to+pv)\ ^ q|_ 

2. ai = inf„_„/3*^(to), Ri(to) > infu,«/3i„(to) and 


y I Hi {to + pu) - Hi {to + pv)\ 
p™ pinfu,^/3L(*o) 


Fi{u,v) 


with {u,v) bounded on [a, 6]^ for some Pi > 0. 

As usually, the proof of weak convergence proceeds in two steps. First, we 
need to show finite dimensional convergence, and then, use a tightness argument. 
Lemma 14.2 and theorem 14.3 in m. for instance, allow then to conclude. 


5.1 Finite dimensional convergence 

As the considered processes are Gaussian, we only have to show the convergence 
of covariance functions. 
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5.1.1 Multifractional Brownian field 

By (O, we compute 


= D [ff(to + pu) + H{^ + pv)] X I|p.(u - y)\\H{to+pu)+H{to+pv) 

+ (2^^(io + pu)\ po + pu\\) X {Hm+ pu) - H{to + pv))^ 

+o (||p.(M - u)|p) +o{H{tQ + pu) - H{to + pv)f (49) 


To show that 


pH(to+pu)+H[to+pv) 

in the neighborhood of p = 0, we study 
[H{tQ + pu)-\-H{tf) + pv) — 2H{to)] Inp = 


p2ff(to) 


H{to + pu) - H(to) 


+ 


\\p.u\r 
H{to + pv) - H{to) 
||p.u||“ 


X ||p.M||“lnp 


X llp.-c^lnp 


for a < Pito). 

As {u\p) I—> ||p.M||“lnp is bounded on [a,b] x [0,1] and 


Vu e [a, 6]; 


H{to + pu) - H{to) p^o 

iip.uih 


we have 


[H {to + pu) + H{to + pv) - 2H (to)] In p 0 
Therefore, in the neighborhood of p = 0, the first term of ijini) is equivalent to 
D[2H{to)]\\u-vf^^*°^ xp2^(*“) 


and the second to 

^ (2iJ(to); ||to||) X {H{to + pu) - H{to + pu))^ 

Let Puv{to) = sup|a;limp^o = o|. We have to distinguish 

the two following cases 

• if H{to) < i-niu,v Puv{to), by definition of /3™(to), 


Vit, u G R 


lim 

p^O 


liJ(to + pu) - H{to + pv)\ 
nHito) 


= 0 


Therefore 


yu,vell^; L;k^(‘«)(p)-y,^(*“)(p)l ^ll[2i/(to)] 

where denotes fractional Brownian field of parameter H{to)- 
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if H(to^ > -jj(to)j 
for all a < in^u,v Puv{to)-, as 


V«,v e Rj; lim Wto + pu) - H(t, + pv)l ^ ^ 

+ P^O 


we have 


Vw, z; e ; -^E [Xt,+pu - Xt,+p^f ^ 0 

Moreover, since there exists u,v G R+ such that H{to) > Puv{to), we can 
consider a G {Puv{to)\ H(to)). The limit 

\H{to +pu) - H{to +pv)\ 

lim sup-= +00 

P^O P°‘ 

implies 

3m, V G R+; lim sup -^E [Xt^+pu - Xt^^+p^f = +oo 

p^O 

Therefore E \Y^{p) — admits a limit for all u,v G R+ when p —> 0 

if and only if 

— inf^z,?; Puviio') 
and 

i: \H(ta+pu) — H(ta+pv)\ -n/.. ^ T)* 

limp^O -i (U, V) G K+ 


In that case, we have for all u, v in R; 


N 

'+ > 


E[Y:{p) - F“(p)]^ ^ (2i7(to); Poll) [r(M,M)]^ 

Remark 6 ITe can see easily that 

hence 

inf /3u{to) < inf/3™(to) 

uGU u,v 


(50) 


(51) 


Conversely, assume there exist u,v GU such that /3u{to) < Pv{to), and let 
aG iPuito); l3y(to)). The inequality 

|g(to + pu) - H{to)\ ^ |H(fo + pw) - g(to + pu)| ^ |g(to + pr;) - g(to)| 


implies 


\H{to +pu) - H{tQ +pv)\ 

hm sup-= +00 

P^O P°‘ 


and therefore a > Puvito)- Then inf„^„/3„^(to) < Pu(to), which 

gives 

iniPuvito) = inf /3„(to) (52) 

u,v uGU 
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5.1.2 Multifractional Brownian sheet 

In the non-isotropic case, using remark |3 consider N independent mBm 
with parameter function Hi 


N 


( 


(^ w 



4 ‘V 

4 

4 VJ 


E [^XtQ^to+pu^Xto,to+pv\ — 

i=l 

As in the isotropic case, for all* G {1,, N}, consider 

^ f r \H,{to + pu) - H,{to + pv)\ ) 

/3™(io) = sup |a; lim- — -= 0| 

Each process is locally asymptoticly self-similar, therefore 


E 


- X - 


‘'0 




p^O 


e[y:^y:^] 


where denotes 

• fractional Brownian motion of parameter Ui = Hiito), in the case Hi(to) < 

• the centered Gaussian process such that 

E[Y^^-Y,^f =Ki, [r,(^r,u)]" 

where at = inf„_„/3J,„(to), in the case Hi{to) > inf„^„/3^„(to) and 


,. \H,{to + pu) - H,{to + pv)\ 

hm -r-;- a.- ^ ^ i{u,v] 

p^O pinG,^/3t„(to) ’ ' 


with Ei bounded on [a,b]^. 
Then we conclude 


E 


AX]l> 

‘■0 ’’•0 




nYi^i 


nYi^i 


p^O 


e[y^y: 


where E“ = H^i Y°‘'- 


5.2 Tightness of laws 

The study of weak convergence is well-known for stochastic processes indexed by 
R+. A comprehensive review was made by Billingsley (cf |H]) for a compact set 
of index ([0,1]). In f[TT]i. Karatzas and Shreeve stated the same kind of results 
for the whole R+. The case of R:5|^ can be found in (cni) whose corollary 14.9 
provides 

Proposition 18 Consider a sequence of continuous processes with 

€ R+l on {n,E,P) such that 
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1. there exists a positive constant v such that 


sup E 

n>l 



V 

< 00 


2. for all T > 0 and all s, t in [0,T]^, there exist positive constants a, (3 
and Ct such that 


supE 

n>l 


XM _ xM 


< CT\\t - s\\^+f^ 


Then the probability measures Pn = P- ' on ,b(c(rI))) 

form a tight sequence. 

We verify the conditions of proposition d in the case of mBm, in the fol¬ 
lowing sections. 


5.2.1 Multifractional Bro-wnian field 

By m, there exist positive constants Kt and Lt such that for all u, v in 
[0,T]^ 

p^^E[Y^{p)-Y,^{p)f = E[Xt,+pu-Xt,+p.f 

< KT\\p.{u-v)f^^*°+P^^ 

+Lt \H{to + pu) - iJ(to + pv)\^ 


Therefore, 


E [F“(p) - Y:{p)f < K'^ ^ ^^ \H{t,Ypu) ^H{to+pv)\ 

• In the case i/(to) < inf„,.i, (3uv{to), there exists Mt > 0 such that 


E 




< Mt \\u — V 


\2H{to) 


In the case E[{to) > inf„_.„ Puvito)^ under the assumption 

\H{to + pu)-H{to + pv)\ 
hm-r-f— 3 — 7-3 -= T{u,v) 


with {u,v) 1 -^ bounded on [a, 6 ]^, there exists Mt > 0 such that 


E 




^ 2 


< Mt IIm — f 


|2(/3Aff(to)) 


Since the process is Gaussian, we get an exponent greater than N in the 
usual way. Then we can conclude by proposition dl that the laws of Y°‘ are 
tight. 
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5.2.2 Multifractional Brownian sheet 


In the same way as in paragraph 14.1.21 

N 


E 


— Y^P'> 


1 


P- 


,2 5:^ a. 


■E 


.i=l 




(^) 


N 






n 


< kY^e 




_ 

po^i 


then, under the assumptions of theorem 0] there exists a positive constant Mt 
such that 


E 


y(p) _ y(p) 


< Mt IIm — 'cll 


We conclude as in the isotropic case. 
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